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Abstract: We discuss the effective theory appropriate for studying soft leptogenesis at 
temperatures T ^ 10^ GeV. In this regime, the main source of the B — L asymmetry is 
the CP asymmetry of a new anomalous i?-charge that couples to generalized anomalous 
electroweak processes. Baryogenesis thus occurs mainly through i?-genesis, and with an ef- 
ficiency that can be up to two orders of magnitude larger than in usual estimates. Contrary 
to common belief, a sizeable baryon asymmetry is generated also when thermal corrections 
to the CP asymmetries in sneutrino decays are neglected which, in soft leptogenesis, im- 
plies vanishing lepton- flavour CP asymmetries. We present general Boltzmann equations 
for soft leptogenesis that are valid in all temperature regimes. 
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1. Introduction 

In the hot and fast expanding Universe, during the first instants after the Big Bang, at 
any given temperature T all particle physics processes having a characteristic time scale 
r larger than the Universe age tu{T) do not occur, and must be neglected. This is im- 
portant, because generically speaking several particle interactions that are allowed by the 
fundamental gauge symmetries violate some other global conservation laws. However, un- 
til the Universe is old enough that these reactions can occur with rates comparable, or 
larger, than the Universe expansion rate, the will-be violated quantities remain effectively 
conserved. In the language of field theory Lagrangians, this means that at each cosmolog- 
ical temperature T, the relevant particle physics processes are determined by an effective 
Lagrangian in which all the parameters responsible for 'slow' reactions, that is reactions 
with characteristic timescales r ^> tu{T), must be set to zero. By doing this, it is then 
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easy to identify the new global symmetries of the effective Lagrangian, and if no anomalies 
are involved, these symmetries correspond to conserved quantities. 

In the context of the early Universe the meaning of 'effective theory', the one that we 
will use in this paper, differs somewhat from what is generally meant in particle physics by 
'effective field theory'. The latter case refers to the low energy theory obtained e.g. from 
a fundamental Lagrangian when all the states heavier than some high energy cutoff are 
integrated out, and corresponds to a theory with a reduced number of degrees of freedom. 
In contrast, the effective theories required to study particle physics processes in the early 
Universe correspond to theories with a reduced number of fundamental parameters. Let 
us explain this in some detail: at each specific temperature T, particle reactions must be 
treated in a different way depending if their characteristic time scale r (given by inverse of 
their their thermally averaged rates) is 

(i) much shorter than the age of the Universe: r ^ tu{T); 

(ii) much larger than the age of the Universe: r » tu{T); 

(iii) comparable with the Universe age: r ~ tuiT). 

The first type of reactions (i) occur very frequently during one expansion time 1/H{T) 
{H(T) being the Hubble parameter at T) and their effects can be simply 'resummed' by 
imposing on the thermodynamic system the chemical equilibrium condition appropriate for 
each specific reaction, that is J2i l^i — YIfI^f-, where /x/ denote the chemical potential of 
an initial state particle, and jjlf that of a final state particle. The numerical values of the 
parameters that are responsible for these reactions only determine the precise temperature 
T when chemical equilibrium is attained and the resummation of all effects into chemical 
equilibrium conditions holds but, apart from this, have no other relevance, and do not ap- 
pear explicitly in the effective formulation of the problem. Reactions of the second type (ii) 
cannot have any effect on the system, since they basically do not occur. Then all physical 
processes are blind to the corresponding parameters, that can be set to zero in the effective 
Lagrangian. In most cases (but not in all cases) this results in exact global symmetries 
that correspond to conservation laws for the corresponding charges, that must be respected 
by the equations describing the dynamics of the system. Reactions of the third type (iii) 
in general violate some symmetries, and thus spoil the corresponding conservation condi- 
tions, but are not fast enough to enforce chemical equilibrium conditions. Only reactions 
of this type appear explicitly in the formulation of the problem (they generally enter into 
a set of Boltzmann equations for the evolution of the system) and only the corresponding 
parameters represent fundamental quantities in the specific effective theory. 

Several examples of the importance of using the appropriate early Universe effective 
theory can be found in leptogenesis studies. Leptogenesis [1,2] was first formulated in the 
so called 'one flavour approximation' [3-5] in which a single SU{2) lepton doublet of an 
unspecified flavour is assumed to couple to the lightest singlet seesaw neutrino, and it is 
thus responsible for the generation of the lepton asymmetry. Until the works in rcfs. [6,7], 
most leptogenesis studies were carried out within this framework, although a few earlier 
works had already explored in some detail the effects of lepton flavours in leptogenesis [8,9], 
or had used them in specific leptogenesis realizations [10]. 
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Nowadays, it is well understood that the 'one flavour approximation' gives a rather 
rough and often unreliable description of leptogenesis dynamics in the regime when flavour 
effects are important. This is because such an 'approximation' has no control over the 
effects that are neglected, and thus the related uncertainty cannot be estimated. On the 
other hand, it is seldom recognized that if leptogenesis occurs above T ~ 10^^ GeV, when 
all the charged leptons Yukawa interactions have characteristic time scales much larger 
than the 'one flavour approximation' is not at all an approximation. Rather, it is 
the correct high temperature effective theory that must be used to compute the baryon 
asymmetry. The corresponding effective Lagrangian is obtained by setting to zero, in first 
place, all the charged lepton Yukawa couplings, so that the only remaining flavour structure 
is determined by the Yiikawa couplings of the heavy Majorana neutrinos, that must remain 
non-vanishing in order that decays into light leptons can occur. 

The transition to the regime where the unflavoured effective theory must be used, 
corresponds to a different physics framework that is characterized by an overall reduced 
amount of CP violation, that is encoded in a single CP violating parameter e, rather than in 
the three (or two, for intermediate temperature regimes [6-8]) flavoured CP asymmetries 
5 e^, €r. Similarly, the lepton density asymmetry is produced into a single species of 
leptons rather than in the three (or two) of the flavoured regimes. Thus, in the regime 
in which the appropriate theory is flavour blind, a lesser amount of baryon asymmetry 
can be produced. Of course, in the same regime also the rates of other processes, like for 
example those induced by the Yukawa couplings of the light quarks [11] or the electroweak 
sphaleron rates [11,12], must be set to zero, but being these typical 'spectator' processes, 
the effect of switching them off is numerically much less relevant.* 

In summary, early works on Standard Model (SM) leptogenesis were carried out from 
the start within the unflavoured effective theory. Quite likely this happened because the 
corresponding Lagrangian is much more simple than the full SM Lagrangian given that 
the number of relevant parameters is reduced to a few. The main virtue of subsequent 
studies on lepton flavour effects [6-9] was that of recognizing that below T ~ 10^^ Qgy the 
unflavoured theory breaks down, and the new theory, that at each step the temperature is 
decreased brings in new fundamental parameters, can give genuinely different answers for 
the amount of baryon asymmetry that is generated. 

In supersymmetric leptogenesis the opposite happened, because the effective theory 
that was generally used is in fact only appropriate for temperatures much lower than the 
typical temperatures T S> 10^ GeV in which leptogenesis can be successful, and only quite 
recently it was clarified that in the relevant temperature range a completely different effec- 
tive theory holds instead [13]. More specifically, it was always assumed (often implicitly) 
that lepton-slepton reactions like e.g. U -f^ U that are induced by soft super symmetry- 
breaking gaugino masses, are in thermal equilibrium (see refs. [2,14,15] for examples of well 
known papers adopting this assumption). This implies equilibration between the leptons 
and sleptons density asymmetries (superequilibration) while in general, in supersymmet- 
ric leptogenesis, superequilibration (SE) does not occur. In fact, requiring that the rates 

* Spectator processes are fast processes that do not violate B — L but that can still have an impact on 
the baryon asymmetry yield of leptogenesis [11, 12]. 
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induced by supersymmetry-breaking scale {Asusy) parameters, like soft breaking masses m 
or the Higgsino mixing parameter n/j, are slower than the Universe expansion rate when 
T ^ M (being M the heavy neutrino mass, and mp below the Planck mass) one obtains 

The effective theory appropriate for studying supersymmetric leptogenesis, in which the 
heavy Majorana masses certainly satisfy the bound Eq. ( |1 . 1| ) , is thus obtained by setting 
m, iifj — )• 0. The consequences of this were analyzed in [13] and are far reaching. At 
T > W GeV, besides the occurrence of non-super equilibration (NSE) effects, additional 
anomalous global symmetries that involve both SU{2) and SU{3) fermion representations 
emerge [16]. As a consequence, the electroweak (EW) and QCD sphaleron equilibrium 
conditions are modified with respect to the usual ones and, among other things, this also 
yields a different pattern of sphaleron induced lepton-flavour mixing [6-8]. In addition, a 
new anomaly-free /^-symmetry can be defined and the corresponding charge, being exactly 
conserved, provides a constraint on the particles density asymmetries that is not present in 
the SM. However, in [13] it was also concluded that, in spite of all these modifications, the 
resulting baryon asymmetry would not differ much from what was obtained in the usual 
scenario. Basically, this happens because by dropping the SE assumption and accounting 
for all the new effects, only modifies spectator processes, while the overall amount of CP 
asymmetry that drives leptogenesis remains the same. 

The most interesting scenario in which the appropriate effective theory not only yields 
far reaching qualitative differences but also very large quantitative effects, is in soft lepto- 
genesis (that is leptogenesis where the origin of CP violation is in the soft supersymmetry- 
breaking terms [17-19]) if it occurs above the SE threshold Eq. ( |1 . 1\ ) . This is because of 
two main reasons: 

(I) The first one is that in soft leptogenesis there is a strong cancellation between 
the CP asymmetries for sneutrino decays into scalars and into fermions e = Es + ej ~ 0. 
This cancellation is almost exact in the T = limit, and e gets lifted to an appreciable 
level only when thermal corrections are included [18, 19]. In the NSE regime however, 
the independent evolution of the scalar and leptonic density asymmetries implies that 
the corresponding efficiencies ij^j are different. When these different 'weights' are taken 
into account, the cancellation between the scalar and fermion contributions to the baryon 
asymmetry gets spoiled and a non- vanishing result is obtained even in the e = Cg + ef ^ 
limit. Note that this effect can dominate over the ones due to thermal (or higher order) 
corrections to the CP asymmetries. This situation is reminiscent of the so called Purely 
Flavoured Leptogenesis (PEL) scenarios [7,20,21,23] where the vanishing of the total CP 
asymmetry resulting from the sum over lepton flavours e = Yla = (q = e, fi, r) does 
not imply a vanishing baryon asymmetry AB oc iJaGa 7^ 0. This, provided that lepton 
flavour equilibrating (LEE) reactions ia ^ that in the PEL case play the same role 
than SE for soft leptogenesis, remain out of equilibrium [22]. 

(II) The second reason is even more interesting. In the high temperature effective 
theory two new global symmetries (a i?-symmetry and a PQ-like symmetry) arise. While 
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these symmetries are anomalous, two new anomaly free combinations of charges involving 
R and PQ can be defined. These new charges, that we denote as Rb and -R^, are only 
(slowly) violated by sneutrino dynamics, that is by reactions of the third type (iii) in 
the classification given above, and thus their evolution must be followed by means of two 
new Boltzmann Equations (BE). Because charge density asymmetries get mixed by EW 
sphalerons, these equations are coupled to the BE that control the evolution of the B — L 
asymmetry, and thus the dynamical evolution oi Rb and R^ affects its final value. What 
is important, is that the CP violating sources for these two charges, that are respectively 
and — e/, are not suppressed by any kind of cancellation, and the corresponding 
density asymmetries remain large during leptogenesis. They act as source terms for B — L 
that is thus driven to comparably large values. As regards the final values Rb and Ry^ 
at the end of leptogenesis, they are instead irrelevant for the computation of the baryon 
asymmetry since, well before the temperature when the EW sphalerons are switched off, 
soft super symmetry-breaking effects attain in-equilibrium rates, implying that R and PQ 
cease to be good symmetries also at the perturbative level. Thus, they decouple from the 
sphaleron processes that then reduce to the usual SM B — L conserving form, that involves 
only quarks and leptons. The baryon asymmetry is then given only hy B — L conversion, 
according to the usual equation B = ^ [B — L). 



The outline of the paper is as follows: In Section ^ we recall the main motivations for 
soft leptogenesis and review recent results and the relevant literature. The soft leptogenesis 
scenario is summarized in Section ^ we first present the relevant Lagrangian and next, 
in Section we compute the various CP asymmetries including some subleading terms 
for the CP asymmetries "in mixing" - generated by sneutrino self-energy diagrams - that 
avoid a complete cancellation between the fermions and bosons contributions even in the 
T — )• limit. Conversely, for the CP asymmetries "in decays" - that are generated by 
vertex corrections ~ we present in Section 3.2 a simple proof ensuring that at one loop 
the cancellation at T = is exact. The effective theory appropriate for studying the 
generation of the baryon asymmetry when the heavy sneutrino masses satisfy the bound 
in Eq. ( |1.1| ) is described in Section We derive the equilibrium conditions and the relevant 
conservation laws that constrain the particle density asymmetries, we identify the new 
quasi-conserved charges, and we also compute the matrices that control the sphaleron 
induced lepton flavour mixing for two different sets of values of the electron and down- 
quark Yukawa couplings. In Section ^ we present the set of the five basic BE, that is valid 
for numerical studies of soft leptogenesis at all temperatures, and in Section 5.1 we discuss 
a simple case in which the role played by the Rb and Ry. charge asymmetries is particularly 
transparent. In Section ^ we compute numerically the amount of baryon asymmetry that 
can be generated in soft leptogenesis within the NSE regime, and compare it to previous 
results based on the assumption of SE. Finally in Section we present a simple explanation 
of the large numerical enhancements, we recap the main results and draw the conclusions. 
Two Appendix complete the paper: in Appendix^ we collect the relevant thermal factors, 
in Appendix p| we present a more complete set of BE which also include scatterings. 
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2. Soft leptogenesis: review and motivations 

With the discovery of neutrino oscillations, leptogenesis [1,2] became a particularly well 
motivated mechanism to explain the cosmic baryon asymmetry. This happened because 
the scale of the oscillation mass square differences is perfectly compatible with sufficient 
deviations from thermal equilibrium in the decays of the heavy seesaw neutrinos [24]. 
For a hierarchical spectrum of the heavy Majorana states, successful leptogenesis requires 
generically a quite large leptogenesis scale [25] , corresponding to seesaw neutrino masses of 
order M > 2.4(0.4) x 10^ GeV for vanishing (thermal) initial neutrino densities [15,25,26].^ 

The presence in the theory of such a large mass scale poses a serious fine tuning problem 
for keeping the Higgs mass parameter at the clcctrowcak scale [31]. Low-energy supersym- 
metry can be invoked to naturally stabilize the hierarchy between this new scale and the 
electroweak one. This, however, introduces a certain conflict between the gravitino bound 
on the reheat temperature and the thermal production of the heavy singlets neutrinos [32] . 

Once supersymmetry is introduced, there are, however, new sources of lepton num- 
ber and CP violation that are related to the soft supersymmetry-brcaking terms involving 
the sneutrinos. This allows for a different leptogenesis scenario that is specific to super- 
symmetry known as 'soft leptogenesis' [17-19]. Given that the new effects are generically 
suppressed by powers of the ratio between the soft supersymmetry-breaking scale and the 
sneutrino masses Agusy/M, the characteristic temperature window in which the new con- 
tributions can give relevant effects is roughly 10^ — 10^ GcV. Thus, apart from a small 
temperature interval lying approximately within 10^ GeV ^ T ^ 10^ GeV, supersymmetric 
leptogenesis can proceed at any temperature above the EW scale, and the low temperature 
soft leptogenesis realization allows to evade the gravitino problem. Note however, that the 
presence of a forbidden temperature window implies that the two different leptogenesis re- 
alizations never overlap. Thus, when studying soft leptogenesis, the standard contributions 
to the CP asymmetries give irrelevant effects and can be safely neglected, and the opposite 
is true when the standard high temperature scenario is assumed. 

In the original papers on soft leptogenesis [18, 19] only one type of contributions to 
the CP asymmetries in sneutrino decays was identified: the so called CP violation in mix- 
ing. CP violation in mixing is induced by the bilinear sneutrino B term that removes 
the mass degeneracy between the two real sneutrino states. As for the case of resonant 
leptogenesis [29], the sneutrino self-energy contributions to the CP asymmetries can then 
be resonantly enhanced, and give rise to sufficiently large CP violation in sneutrino de- 
cays. However, to satisfy the resonant condition, unconventionally small values of B are 
required [18,19,33,34]. Extended scenarios were thus proposed in order to alleviate this 
problem [35,36]. However, it was also realized that within the context of the minimal 
scenario, besides CP violation in sneutrino mixing additional sources of CP violation can 
arise from vertex corrections to the decay amplitudes, and from the interference between 
mixing and decay [37,38]. These new sources of CP violation (the so called "new ways to 
soft leptogenesis" [37]) are induced by gaugino soft supersymmetry-breaking masses that 

t These limits are basically not affected by flavour effects [27,28]. With resonantly enhanced CP asym- 
metries [29] or in various extended scenarios [30] lower leptogenesis scales are instead possible. 
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appear in vertex corrections to the decays. With respect to the mixing contributions, these 
corrections are suppressed by more powers of Asusy/M, and thus they can be sizable only 
at relatively low temperatures T ^ 10^ GeV [37] . Although in this regime they can allow 
for more conventional values of B, such a low seesaw scale implies that the suppression 
of the light neutrino masses is mainly due to very small values of the Yukawa couplings 
~ lO""*^" rather than to the seesaw scale. 

Concerning the role of flavour [6-9,22,27,28,39-41] and spectator effects [11,12], they 
had all been neglected in the original soft leptogenesis papers [18,19,37] that were based on 
the single-flavour scenario. However, soft leptogenesis always occurs at temperatures where 
the appropriate effective theory must include the effects of the three lepton flavours. This 
was done in Ref. [33] that also included spectator effects, but that assumed a constrained 
scenario with universal trilinear couplings. Within this context, it was found that the 
leptogenesis efficiency could be enhanced by 0(30) with respect to the single flavour anal- 
ysis. The more general scenario of non-universal trilinear couplings was considered in [42], 
that also included the possibility of damping flavour effects through large LFE spectator 
processes [22]. It was found that when the assumption of universality is dropped, flavour 
effects can play an even more important role, with the possibility of enhancing the lepto- 
genesis efficiency by more than three orders of magnitude with respect to the one flavour 
treatment. 

In spite of all these advancements and refinements in soft leptogenesis studies, a cru- 
cial point has been always overlooked. As was first pointed out in [13], when the sneutrino 



masses satisfy the bound Eq. (LI) the appropriate effective theory for studying early Uni- 
verse processes in a supersymmetric scenario is different from the one that has always been 
used. As we will show, in soft leptogenesis the correct effective theory implies particularly 
dramatic effects, namely the final baryon asymmetry that is produced can be up to two 
orders of magnitude larger than what is obtained with the previous treatments. 

3. Soft Leptogenesis Lagrangian and CP asymmetries 

The superpotential for the supersymmetric seesaw model is: 



W = -MijN^N^ + Yi^N^l^H^, (3.1) 



where i,j = 1,2,... label the chiral superfields of the heavy SU{2) singlet Majorana 
neutrinos defined according to usual conventions in terms of their left-handed Weyl spinor 
components {N'^ has scalar component A^* and fermion component Nf^), a = e,fi,T labels 
the flavour of the SU{2) lepton doublets £a = {va,&a)^ ^ = {H^^H^Y^ denotes the 
up-type Higgs doublet superfield, and contraction of the SU (2) indexes between doublets 
laHu = ^pa^aHu with £12 = +1 is left understood. 

The relevant soft supersymmetry-breaking terms involving the scalar components of 
the A^*^ superfields and the SU (2) gauginos A2 are given by 

-Csoft = MfjN*Nj + (^AYi^N*£^H^ + ^BMijN^N* + ^tuXPl^ + h-c j . (3.2) 
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C/(l) gauginos can be straightforwardly included in similar form. In Eq. (|3.2| ) we have 
assumed for simplicity universal trilinear and bilinear couplings Aia = AYia and Bij = 
BMij. The sneutrino and anti-sneutrino states mix, resulting in the mass eigenstates: 

v2 

N_i = ^(e**/27V* - e-'*/2^,), (3-3) 
V2 

where $ = arg(i?M), and N±i have mass eigenvalues: 

Ml, = Mf, + Mf,±\BMu\. (3.4) 

The interaction Lagrangian involving the mass eigenstate sneutrinos N±i, the Ma- 
jorana singlet neutrinos Ni, the SU{2) gauginos A2 and the (s)leptons and Higgs(inos) 
doublets reads: 

N+i (HuPL^a + {A + Mi)I^H^) + iN-i (H^PLia + {A- Mi)IaHu 



Y 



+ Yi 



la 



H-PlNda + NiPLiaHu 



+92 



xtPiia CT±C + xtpRKa^Hu - {x^PL^o^T^t^ + 'x^PrKg-.H^ 



+h.c., (3.5) 



where Pl,r are respectively the left and right chiral projectors, (t± = (cji it icj2)/2 with cxj 
being the Pauli matrices, and SU{2) contractions like cr± £^ = (c±)pcr^a* are again 
left understood. All the parameters appearing in the superpotential Eq. ( |3.1[ ) and in the 
Lagrangian Eq. (^]^) (and equivalently in the first two lines of Eq. ( |3.5D ) are in prin- 
ciple complex quantities. However, superfield phase redefinition allows to remove sev- 
eral complex phases. Here for simplicity, we will concentrate on soft leptogenesis arising 
from a single sneutrino generation i = 1 and in what follows we will drop that index 
(Yq, = Yia, N± = N±i, etc.). After superfield phase rotations, the relevant Lagrangian 
terms restricted to i = 1 are characterized by only two independent physical phases: 

cj)A = arg{AB*), (3.6) 

(1^9 = ^arg(Sm^), (3.7) 

which we choose to assign respectively to the slepton-Higgs-sneutrino trilinear soft break- 
ing terms, and to the gaugino coupling operators respectively. In what follows we will 
keep track of these physical phases explicitly and, differently from the convention used 
in Eqs. (|3.l|), (|3.2[) and (|3.5D , we will leave understood (unless when explicitly stated in 
the text) that all the other parameters Ya,B, m2, A etc. correspond to real and positive 
values. 

3.1 CP asymmetries in soft leptogenesis 

Neglecting super symmetry-breaking effects in the heavy sneutrino masses and in the vertex, 
the total singlet sneutrino decay width is given by 

-^N. = ^- ^Z^^- = ' (3.8) 
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where Vu = v sin/3 (with i;=174 GeV) is the vacuum expectation value of the up-type Higgs 
doublet and meg = Yla ^a'^t/^ is the rescaled decay width, that is related to the washout 
parameter K as K = T^/H{M) = mes/'m'*i where the equilibrium mass is defined as 

~ \/ ^ ^Tn^ ~ 10"'^ eV with g* the total number of relativistic degrees of freedom 
(g* = 228.75 in the MSSM). 

There are three contributions to the CP asymmetry in N± decays into fermions 
{Hu, la) and other three for decays into bosons (i/^, £a) [37, 38]. Denoting the total 
CP asymmetries into the scalar and fermion channels respectively with s and / subscripts, 
they are: ef j arising from self-energy diagrams induced by the bilinear B term; e^j arising 
from vertex diagrams induced by the gaugino masses; j arising from the interference of 
self-energy and vertex diagrams. They can be written as: 



ef(r) _ A 4sr 



A/(r) M4B2 + r2 

ef(T) _ A ABT 




sin(/)A, (3.9) 

(3.10) 



(T) eY (T) 3a2 A 1712 rn^ \ ( , ^ , n ^ . m . ^ 

= = ^Mli ['^^^f^^) {-(^- + 20,) - -sm(20,) } (3.11) 



4 (T) _ ei (T) _ 3a2 A ma 



A/(r) A, (T) 2 M M AB^ + \ mj + M'^ 



2 , sin(/)^cos(2(/)g) , (3.12) 



-2 



where a2 = To take into account the effects of lepton flavours, we need to define instead 
of Eqs. ( |3.9D -( ^T2|) the corresponding asymmetries for N± decays into fermions and scalars 
of a specific lepton fiavour £q, and la- The assumption of universality of the soft terms 
implies that the flavour CP asymmetries are simply related to the total CP asymmetries 
through the corresponding decay branching fractions, denoted by Pa'- 

S,V,I _p S,V,I .O^gX 

where, in terms of the Yukawa couplings, the branching fractions can be written as: 

E^" = l- (3-14) 

/3 

It is worth recalling that when the condition of universality of the soft terms is dropped 
the simple relation Eq. ( |3.13| ) does not hold anymore. The expressions for the flavoured 
CP asymmetries in the general case of non-universal soft terms can be found in [42]. The 
Asj(r) terms in Eqs. ( |3.9| )-( |3rT^ ) denote the scalar and fermion thermal factors that are 
related to thermal phase-space, Bose enhancement and Fermi blocking. They are the same 
for N± and are normalized so that their zero temperature limit is Asj{T = 0) = i. Their 
explicit expression is given in Appendix 0. Note that in writing Eqs. (|!|)-(|1|) we have 
implicitly assumed that the thermal factors are flavour independent. This is indeed an 
excellent approximation as long as zero temperature slepton masses and small Yukawa 
couplings are neglected. 
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Eqs. ( |3.9| )- (|3.12| ) are approximate expressions with only the leading terms in 



M 

included. In Eqs. (|3.9|) and (|3.10|) however, we have included also O ( ^^a" ) corrections, 
since when the resonant condition F ~ 2B is satisfied e*^ is the dominant CP asymmetry. 
Note that with these corrections included the self-energy CP asymmetries for fermions and 
scalars do not cancel anymore in the T — )• limit. In contrast, as we will argue in the 
next section, for the vertex and interference asymmetries e^'^ at one loop the cancellation 
between scalar and fermion contributions is exact. Neglecting for simplicity the higher 
order terms in Eqs. ( |3.9| ) and ( |3.1C1| ), the total CP asymmetry summed over scalars and 
fermions final states can be written as 

ea{T) = Pae-[AsiT)-Af{T)] (3.15) 

where e is independent of T and 

[A, -A;] ^-^ 0. (3.16) 
3.2 Vanishing of the CP asymmetry in decay 

The new sources of direct CP violation from vertex corrections involving the gauginos were 
first introduced in soft leptogenesis in [37]. In the same paper it was also stated that the 
new contributions do not require thermal effects to produce a sizable lepton asymmetry in 
the plasma. Gaugino contributions to soft leptogenesis were reconsidered in [38] where it 
was instead found that the zero temperature cancellation between the CP asymmetries for 
decays into scalars and into fermions holds also when vertex corrections are included. This 
issue is of some interest, because if thermal corrections are necessary for soft leptogenesis 
to work, then non-thermal scenarios, like the ones in which sneutrinos are produced by 
inflaton decays and the thermal bath remains at a temperature T <^ M during the following 
leptogenesis epoch, would be completely excluded. Here we present a simple but general 
argument proving that the direct leptonic CP violation in sneutrinos decays vanishes at 
one loop, due to an exact cancellation between the scalar and fermion contributions, in 
agreement with the explicit calculation in [38]. We should also clarify that while above the 
limit in Eq. ( [1.1| ) soft leptogenesis can be successful even when the T — t- limit for the 
decay CP asymmetries is taken, this happens because of other effects which are linked to 
the washout processes, and thus do require a thermal bath. Therefore, the fact that soft 
leptogenesis cannot work in non-thermal scenarios is always true. 

Let us take for simplicity $ = in Eq. (|3.3|) (this amounts to assign the phases (f>A and 
(f)g in Eq. (| 



,7V (TV*) 



and Eq. (3/7) respectively to A and m2), and let us introduce for the various 
amplitudes the shorthand notation Af = A{N± — >• £Hu), A 
with similar expressions for the other final states. From Eq. 

2 



A{N{N*) IHu 



I) we can write 





2 




2 

+ 






2 




2 








4^ 


+ 


aN* 



±2Re(<-Af) 



±2Re(^f • A 



(3.17) 
(3.18) 

where the complex conjugate amplitudes in the last terms of both these equations have 
been rewritten as follows: {Aff = At = Af and {Af')* = A^~ = Af by using CPT 
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Figure I: Soft leptogenesis diagrams for sneutrino decays into scalars (la), (16), (Ic) and into 
fermions (2a), (26), (2c). 



in the second step. The direct CP asymmetry for Nj- decays into fermions is given by the 
difference between Eq. (3.17) and Eq. ( p.l8| ): 



A 



±|2 



Ai 



A 



N* 



+ 



A 



N 



A 



N 



(3.19) 



With the replacements i ^ i and I ^ I*, a completely equivalent expression holds also for 
the decays into scalars. 

The tree level and one loop diagrams for the various decay amplitudes into scalars and 
fermions are given in Figure |. We note at this point that A~ has no one-loop amplitude 
to interfere with (see diagram (1 a)) and thus, up to one-loop, the full amplitude coincides 
with the tree level result, and is CP conserving. A^ is a pure one- loop amplitude (see 
diagram (2 c)) and therefore is also CP conserving. This implies: 



(3.20) 
(3.21) 





2 






2 




aN 







We can thus change simultaneously the signs of 



A 



N 



and 



A 



N* 



m Eq. (|3l9| ) without 



affecting the equality, and the same we can do in the analogous equation for the scalars. 
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This gives: 



\A 



±|2 



At 



A 



N* 



+ 



+ 



A 



N* 



A 



N 



A 



N* 



A 



N 



+ 



+ 



A 



N 



A 



N 



(3.22) 
(3.23) 



Using CPT ^ffi* = A^i^ and unitarity 



N 



AN 
^all 



we can readily see that the sum of 



all - -]\f 

these two equations vanishes. We have thus proved that for and N- independently, 
at one loop there is an exact cancellation between the scalars and fermions final state 
contributions, and thus at T = the direct decay CP asymmetries vanish. 



4. Soft leptogenesis above the superequilibration temperature 

We now discuss the early Universe effective theory appropriate for studying soft leptogen- 
esis in the regime in which superequilibrating reactions like i£ ^ ii, that are induced by 
gaugino masses, and higgsino mixing transitions, that are induced by the supersymmetric 
fi^HuHd term, do not occur. For simplicity, we assume equal masses for all the gauginos 
nil = = = and that the supersymmetric higgsino mixing term also has approx- 
imately the same value: f^H — f^g = ^susy The regime we are interested in is defined by 
the condition given in Eq. (|l.l| ), that is the lower limit on the relevant temperatures is: 

4.1 Anomalous and non anomalous symmetries 

The supersymmetric effective theory appropriate to study particle physics processes in 
the early Universe when the thermal bath temperature satisfies the condition Eq. ( [4.1| ) 
is obtained by setting mg,^^ — )• [16]. In this limit the theory gains two new U{1) 
symmetries: ;Ujj — )• yields a global symmetry of the Peccei-Quinn (PQ) type, and by 
setting also rUg — )• one additional global i?-symmetry arises. 

The charges of the various states under R and PQ, together with the values of the 
other two global symmetries B and L are given in Table |l|. Like L, also R and PQ are not 
symmetries of the seesaw superpotential terms MN^N^ + XN^iHu, since it is not possible 
to find any charge assignment that would leave both terms invariant. In Table || we have 
fixed the charges of the heavy A^'^ supermultiplets in such a way that sneutrinos do not 
carry any charge. This has the advantage of ensuring that all the sneutrino bilinear terms, 
corresponding to the mass parameters M, M, B, are invariant, and thus sneutrino mixing 
does not break any symmetry. However, since R{N^N^) = 0, it follows that the mass term 
for the heavy Majorana neutrino breaks R by two units. ^ 
'"This differs from the assignments adopted in Ref. [13]. 

^Under i?-symmetry the superspace Grassmann parameter transform as — >■ e^^d . Invariance oi J dOd = 
1 then requires R{d6) = — 1. Then the chiral superspace integral of the superpotential J d9^ W is invariant 
if R{W) — 2. By expanding a chiral supermultiplet in powers of 6 it follows that the supermultiplet R 
charge equals the charge of the bosonic scalar component R{b) = R{f) + 1, and thus for the fermion bilinear 
term ii(iV|iV£) = -2. 
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Q 




d" 


e 




Hd 


Hu 


N" 


B 





1 

3 


1 

3 


1 

3 

















L 














1 


-1 











PQ 








-2 
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-1 
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-1 


2 





R 1 

b 


1 


-1 


-3 


1 


-1 


1 


-1 


3 


-1 


2 





-2 


2 





2 





4 






Table 1: B, L, PQ and R charges for the particle supermultiplets that are labeled in the top 
row by their L- handed fermion component. Note that we use chemical potentials for the R- handed 
SU{2) singlet fields u, d, e that have opposite charges with respect to the ones for u'^, d'^, e*^ given 
in the table. The R charges for bosons are determined by R{b) — R{f) + 1. 



All the four global symmetries B, L, PQ and R have mixed gauge anomalies with 
SU{2), and R and PQ have also mixed gauge anomalies with SU{3). Two linear combi- 
nations of R and PQ, having respectively only SU{2) and SU{3) mixed anomalies, have 
been identified in Ref. [16]. They are: ^ 

R2 = R-2PQ (4.2) 
R3 = R-3PQ. (4.3) 

The values of i?2,3 for the different states are given in Table |^. The authors of Ref. [16] have 
also constructed the effective multi-fermions operators generated by the mixed anomalies: 

Oew = n„ {QQQ£a) HuHd , (4.4) 
Oqcd = n, iQQu'^d^ f . (4.5) 

Given that we have three charges R2, B and L with mixed SU{2) anomalies, it is then 
possible to define two anomaly free combinations. The most convenient are B — L and 

Rb = ^B + R2, (4.6) 

whose values are also given in Table For the problem at hand, i?^ is more convenient 
than the charge TZ = |-B — L + R2 = Rb — {B — L) that was introduced in [13]. This is 
because TZ is not conserved in sneutrino decays that are induced by sneutrino-related soft 
terms, and so it does not correspond any more to a global neutrality condition [13]. On the 
other hand, the fact that Rb does not contain any B — L fragment, ensures that it will not 
enter in the final computation of the baryon asymmetry that will only depend on B — L. 
The fact that Rb is independent of L renders also easier writing a BE for its evolution. 

The Rb values in Table |2| imply that the superpotential term N'^ iH^ has charge 
Rb = 2 and thus is invariant. It follows that sneutrinos decays into fermions conserve 

^With respect to Ref. [16], for definiteness we restrict ourselves to the case of three generations Ng = 3 
and one pair of Higgs doublets Nh = 1, and we also normalize _R2,3 in such a way that R2,3{b) — R2,3{f) + 1- 
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-2 


2 








R, ' 
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-1 
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-2 
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-5 
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-3 


-1 
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-1 
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-4 


3 


-2 





Rb ^ 
b 


1 


7 
9 


7 
9 


11 

9 


1 


-3 


1 


-1 


-1 


2 


2 
9 


16 
9 


2 
9 


2 


-2 


2 









Table 2: Charges for the fermionic and bosonic components of the SUSY multiplets under the 



i?-symmetries defined in Eqs.(4.2), (4.3) and (4.6). Supermultiplets are labeled in the top row by 
their L-handed fermion component. We use chemical potentials for the R-handed SU{2) singlet 
fields u, d, e that have opposite charges with respect to the ones for u'^, (F , e'^ given in the table. 



Rb- In contrast, the soft A term in Eq. (3.2) responsible for sneutrinos decays into scalars 



violates Rb by 2 units, more precisely N± — )• Huf^ has /S.Rb = +2 while N± — )• H^H* has 
/S.Rb = —2. As regards the heavy neutrinos, their mass term violates Rb by two units. 
Note that this is precisely like the case when one chooses to assign a lepton number —1 to 
the singlet neutrinos A^. Accordingly, the decays of the heavy Major ana neutrino violate 
Rb by one unit: A^ — )• IHu, iHu have ARb = +1 and the decays to the CP conjugate 
states have ARb = — 1- Since all Rb violating reactions have, by assumption, rates that 
are comparable to the Universe expansion rate, the evolution of this charge must then be 
tracked by means of a specific BE. 



At temperatures satisfying the condition Eq. (4.1) there is at least one other anomalous 
global symmetry, that we will denote by x- It corresponds to U{1) phase rotations of the 
u"^ chiral multiplet that, for its fermionic component, can be readily identified with chiral 
symmetry for the right-handed up-quark. In fact, above T ~ 2 x 10^ GeV, reactions 
mediated by do not occur and the condition hu ^ must be imposed, resulting in 
a new anomalous 'chiral' symmetry. In the SU (3) sector we then have two anomalous 
symmetries i?3 and X) and one anomaly free combination can be constructed. Assigning to 
the L-handed supermultiplet a chiral charge x = ~1 this combination has the form [13] 

Rx = ^ui + Rs, (4.7) 

where = 1/3. When the additional condition /i^ — )• is imposed, a chiral symmetry 
arises also for the d"^ supermultiplet. A second anomaly free Ry. symmetry can then be 
defined in a way completely analogous to Eq. (|^, with 

Kd", — K'u'i — 1/3 [13]. As regards 
perturbative violations of R^, this charge inherits the same violations R3 suffers. The soft 
A term in Eq. (|3.2D violates R3 by one unit, and so do sneutrinos decays into scalars. 
Moreover, since N'^ iHu has an overall charge R3 = 1, a violation by one unit occurs also 
for sneutrinos decays into fermions. Correspondingly, we have AR^ = +1 for the decays 
N, N* hJ, HJ and AR3 = -1 for N, N* HJ, H*t. Of course, similarly to Rb, 



-14- 



also the evolution of Ry. needs to be tracked by means of a BE. 
4.2 Chemical equilibrium conditions and conservation laws 

Because of the network of fast particle reactions occurring in the thermal bath, asymmetries 
generated in sneutrino decays spread around among the various particle species, and this 
can affect directly or indirectly leptogenesis processes. In principle there is one asymmetry 
for each particle degree of freedom. There are however several conditions and constraints 
that reduce the number of independent asymmetries to a few. The three types of reactions 
that have been classified in the introduction give rise to three different types of constraints 
and conditions, that need to be formulated in their own appropriate way: 

(i) Constraints imposed by reactions whose rates are much faster than the Universe 
expansion have to be formulated in terms of chemical equilibrium conditions for the 
chemical potentials of incoming m and final state particles hf- 

Y,^ii = Y.^'F- (4-8) 

/ F 



(ii) Conservation laws that arise when all the reactions that violate some specific charge 
are much slower than the the Universe expansion have to be formulated in terms of 
particle number densities An = n — n and, for a generic charge Q, read: 

Q = ^ Qi/^rii = const, (4.9) 

i 

where Qi is the charge of the i-particle species. We will always assume as initial 
conditions for leptogenesis that all particle asymmetries vanish, and thus we will put 



the constant value of Eq. (4.9) equal to zero. 



(iii) Reactions with rates comparable with the Universe expansion have to be treated by 
means of appropriate dynamical equations. In this case, in order to reabsorb the 
dilution effects due to the Universe expansion, it is convenient to introduce as basic 
variables the number densities of particles per degree of freedom g normalized to the 
entropy density s: 

Fa, = 4.10 

9i s 

Clearly, /ij, An^ and Y^. are all related to particles asymmetries. In particular, the 
number density asymmetries of particles for which a chemical potential can be defined are 
directly related with this chemical potential. For both bosons (6) and fermions (/) this 
relation acquires a particularly simple form in the relativistic limit <C T, and at first 
order in fJ^bj/T <^ 1: 

An, = |rV6, An^ = ^rV/. (4.11) 

While we will always express the various constraints using the most appropriate quantities, 
eventually to solve for the large set of conditions in a closed form we will need to use a 
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single set of variables. We will take this to be the set {YAi}) and will leave understood 
that our solutions to the constraining conditions are obtained after expressing /ij and Anj 
in terms of this set, through Eq. ( [4.11 ) and Eq. ( 4.10| ). 



In the following we denote the chemical potentials with the same notation that labels 
the corresponding field: (p = fi^ and, for definiteness, we fix the relevant values of the 
temperature around T ~ 10*^ GeV. The set of conditions that constrain the particle abun- 
dances at this temperatures are listed below, more details about the various constraints 
can be found in [13]. 

1. At T » Mw, gauge fields have vanishing chemical potential W = B = g = ^ [43]. 
This also implies that particles belonging to the same SU{2) or SU (3) multiplets 
have the same chemical potential. For example cj){l3 = +i) = (/)(/3 = — ^) for a field 
(j) that is a doublet of weak isospin /, and similarly for color. 

2. Denoting by Wr, Bji and gji the right-handed winos, binos and gluinos chemical 
potentials, and by i, Q {i, Q) the chemical potentials of the (s)lepton and (s)quarks 
left-handed doublets, the following reactions: Q+gR — )■ Q, Q + Wr — )■ Q, I+Wr — >• £, 
I -|- Br — >• ^, imply that all gauginos have the same chemical potential: 

-~g = Q-Q = -W = l-i=-B, (4.12) 

where we have introduced W ^ B and g to denote the chemical potential of the left- 
handed gauginos. It follows that the chemical potentials of the SM particles are 
related to the chemical potential of their respective superpartners as 

Qj=QJ^~g (4.13) 
Hu,d = Hu4+~9 (4.14) 
d,e = u,d,e — g. (4-15) 



u 



The last relation, in which u,d,e = UR,dR,eR denote the i?-handed SU{2) singlets, 
follows e.g. from u'f^ = u'j^ + g for the corresponding L-handed fields, together with 
u'l = —ur, and from the analogous relation for the SU{2) singlet squarks. 



Eqs.(4.13)-(4.15) together with the vanishing of the chemical potentials of the gauge 
fields and the equality of the chemical potentials for all the gauginos, implies that we are 
left with 18 chemical potentials (or number density asymmetries) that we chose to be the 
ones of the fermionic states. They are 15 for the SM quarks and leptons, 2 for the up- 
type and down-type higgsinos, and 1 for the gauginos. These 18 quantities are further 
constrained by additional conditions. 

3. Before EW symmetry breaking hypercharge is an exactly conserved quantity. There- 
fore for the total hypercharge of the Universe we have 

ytot = ^ An^y;, -h ^ An/y/ = 0, (4.16) 
b f 
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where y^j denotes the hypercharge of the 6-bosons or /-fermions. It is useful to 
rewrite exphcitly this condition in terms of the rescaled density asymmetries per 



degree of freedom {Ia^} defined in Eq. ( 4.10 ): 



J] (Yaq. + 2yA«, - lAdJ - E (^A£. + ^"AeJ + Y^f,^ - Y^^^ = 0. (4.17) 

i a 

4. Chemical equilibrium for reactions that are mediated by the leptons and quarks 
Yukawa couplings give: 

^„ -e„ + #rf + 5 = 0, (a = e, /i, r), (4.18) 
Qi-di + Hd+~g = {di = d,s,b), (4.19) 

Qi-Ui + Hu + g = 0, {ui = c, t). (4.20) 

At T ~ 10^ GeV, Yukawa equilibrium for the up quark is never realized. For a = e 
and for the d-quark Yukawa equilibrium holds as long as T ^ 10^(H-tan^ /3) GeV [44] 
and T < 4 • 10^(1 + tan^ /?) GeV respectively. Then, for T ~ 10*^ GeV both condition 
hold only if tan (3 ^ 35, while they both do not hold if tan /? ^ 5. As we will discuss 
below, in the latter case the Yukawa equilibrium conditions get replaced by other two 
conditions, and thus the overall number of constraints does not change. Below we 
present results for the large and small tan /3 cases, and since they do not differ much, 
we omit the corresponding results for the intermediate case 5 ^ tan (3 ^ 35. 
Besides the previous Yukawa equilibrium conditions, quark intergenerational mixing 
guarantees that below T ^ 10^^ GeV the three quark doublets have the same chemical 
potential: 

Q = Qs = Q2 = Qi. (4.21) 

5. Finally, reactions induced by the QCD and EW sphaleron multi-fermion opera- 
tors Eq. (U) and Eq. (|]|) imply [16] 

3Y^Qi + ^ea + Hu + Hd + 4.g = 0, (4.22) 

i a 
i i 

Counting the number of additional conditions listed in items 3 to 5, we have 1 from 
global hypercharge neutrality, 8 from Yukawa equilibrium plus 2 due to quark intergener- 
ational mixing, and 2 from the EW and QCD sphaleron equilibrium. This adds to a total 
of 13 constraints for the initial 18 variables, meaning that 5 quantities must be determined 
from dynamical evolution equations. These quantities can be chosen, for example, as the 
density-asymmetries of the three lepton flavours Y^i^ , of the up-type higgsinos Y^jj and 
of the gauginos Y^g, where the last one allows to relate the previous four quantities to the 
corresponding densities asymmetries of their superpartners. This choice would be a natural 
one since these are the density asymmetries that 'weight' the various interactions entering 
the BE for soft leptogenesis. However, the EW and QCD sphalerons reactions Eq. ( [4.4| ) 
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and Eq. ( ^.51) imply fast changes of these asymmetries. A much more convenient choice 
is instead that of using appropriate hnear combinations of the various asymmetries cor- 
responding to anomaly free and quasi-conserved charges, where with 'quasi-conserved' we 
refer to charges that are not conserved only by the 'slow' sneutrino-related reactions. These 
quantities can be identified with the three flavoured leptonic charges Aq, = B /3 — and 
with the two Rb and charges discussed in the previous section. In terms of the rescaled 
density asymmetries per degree of freedom they read: 

Ya^ =6Yaq+Y^ {Yau. + YAd, ) - 3 {2YAi,, + ^Ae, )-2YAg, (4.24) 

i 

YaRs = -6Yaq - (131^A«, -5YAdJ 

i 

+ 5^ (10 1a£. + 7 Yas^ ) + 68 YAg + 10 Y^^^^ - 2 F^^^ , (4.25) 



Yar^ = 3 (3 - 2 YAg) + - Yar, , (4.26) 



1 

3 

where, in the last expression. 



Yar, = -ISYaq - 3 J] (11 Yau, -^YAdJ 

i 

+ E + 13lAeJ + 82YAg + IQY^fj^ - UY^j,^. (4.27) 



The density asymmetries of the five charges in Eqs. (|4.24D -( ^!26|) then define the basis 



Ya^ = {Yao, ) YaRb , Yar^ } in terms of which the five fermionic density-asymmetries Ya^p^ = 
{YAia^ YAg, Y^fj^}, that are the relevant ones for the soft leptogenesis processes, have to 
be expressed. We will do this by introducing a 5 x 5 A-matrix defined according to: 

Ya^^ = Aab Ya, , (4.28) 



where the numerical values of Aab Eire obtained from Eqs. ( [4.24D -( p:.26 ) subjected to the 



constraining conditions listed in items 3 to 5. Let us note at this point that the 3 x 
5 submatrix vl^^b for the lepton densities represents the generalization of the A matrix 
introduced in [8], A^^^^ generalizes the Higgs C-vector first introduced in [11], and Agi, 
generalizes the C-vector for the gauginos first introduced in [13]. As regards the density 
asymmetries for the bosonic partners of ia and of H^, they are simply given by: A^ ^ = 

2 {Ae^ b + Agb) and 6 = 2 (^Ag^ + Agb^ . 

4.3 Case I: Electron and down-quark Yukawa reactions in equilibrium 

If the down-type Higgs vev is relatively small v^, ^ v, the values of the electron and 
down-quark masses are obtained for correspondingly large values of the and hg Yukawa 
couplings. For v^/vd = tan /3 ^ 35 we have a regime in which at T ~ 10*^ GeV, that is 
well above the NSE threshold Eq. ([4.1| ), both and related reactions are in equilib- 
rium. In this case all the eight Yukawa conditions Eqs. ( 4.18D -( p:.20 ) hold. Solving for 
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the densities-asymmetries Y^^^ = {Ia^q) ^Ag, ^A_ff„} terms of the charge-asymmetries 
= {^Aa5^AKB>^Ai?x} subject to the constraints in items 3 to 5, yields 
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12798 
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-102411 


-102411 


-102411 


108108 


-335907 / 



(4.29) 



4.4 Case II: Electron and down-quark Yukawa reactions out of equilibrium 

If Vd is not much smaller than Vu, resulting in tan/3 ^ 5, then both and are sufficiently 
small that at T 10^ GeV the related Yukawa reactions do not occur. In this case we have 
to set hd, /le — )• and the corresponding two Yukawa equilibrium conditions in Eqs. ( [4.18 )- 
( |4.19 ) do not hold. However, two conservation laws replace these conditions, /ig — 
implies that we gain a 'chiral' symmetry for the right-handed fermion and scalar electrons, 
ensuring that the total number-density asymmetry Ang + Ang is conserved. As usual, we 
assume that the constant value of this quantity vanishes, which in terms of the rescaled 
density asymmetries per degree of freedom implies: 



>Ao = . 



(4.30) 



For the right-handed down quark we could define an anomaly-free charge completely equiv- 
alent to YAiJx in Eq. (|^ but, given that in this regime all the dynamical equations are 
symmetric under the exchange u -f-T- d, it is equivalent, and much more simple, to impose 
the condition 

Y^d = Y^u . (4.31) 

The net result is that, with respect to the previous case, the total number of constraints is 
not changed, and again five quantities suffice to express the rescaled density asymmetries 
for all the fields. For the 5x5^ matrix defined in Eq. ( 4.28| ) we obtain: 



A 



1 



9 X 162332 
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-210531 21573 
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8678 -3197 


-165529 
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7497 7299 


7299 


23634 


4833 
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-11322 -18477 


-18477 


23940 


-74385 / 



(4.32) 



5. Basic Boltzmann Equations 

In order to render clear the role played by the new charges ARb and AR^ and by NSE 
effects, in this section we introduce a simplified set of BE including only decays and inverse 
decays of heavy neutrinos and sneutrinos. However, for the numerical results that are dis- 
cussed in the next section, we have used the more complete (and involved) set of equations 
described in Appendix |B|. 
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The evolution of the number density of the heavy states normahzed to the entropy 
density s is given by 



Yn-- 




(5.1) 
(5.2) 



where the time derivative is defined as Y 



sHz^ with z 



M/T, and H = H{z) is the 
Hubble parameter. In Eq. (|5.1| ) 7Ar represents the (thermal averaged) total decay width of 



the heavy neutrino into particles and sparticles of all a-flavours 7Ar = 7^ and Y^ 
the equilibrium density. For the heavy sneutrinos, we denote with N the sum of 



and A^_. Thus in Eq. (5^) = +^jv ' '^'^i^^ Y~ represents the equilibrium density 



"Ar+ ' "AT-' 

of a single sneutrino. For the reaction rates we have 



(5.3) 



p=sj a 



where the p sum in the r.h.s of the last equality is over s-scalars and /-fermions final states, 
while 7^_^ = 7^_ = 7^/2. 

In writing down the evolution equations for the five charges Ya^ ■, YaRb j Yar^ it is 
convenient to introduce a special notation for the scalars and fermions density asymmetries 
(per degree of freedom) normalized to the respective equilibrium densities Yg'^ = 2Yf'^ = 

15 . 



yAs,Af 



YAs,Af 



(5.4) 



Using Eqs. ( 1.11 ) and ( 4.10 ) together with ( 4.13 ) and ( 4.14 ) it is then easy to verify that 



y 



y. 



M,AHu 



+ yA. 



9 ■ 



(5.5) 



Including only decays and inverse decays, the Boltzmann equation for the flavour charges 
Aq, = B /3 — La read: 



Yac 




2]^ + [yM^+y^f^. 



2 1 ^ + {y^^ +yAH. 



7^ / X 

+ [yAi^ + yAHu) 



7 



N 



+ {y^i+y. 



AH, 



. (5.6) 



To write this expression in a more compact form, we define the total flavoured CP asym- 
metry e" = + and the total sneutrinos decay rate into a leptons and sleptons 



'n 



7~^" + 7^°- For quantities without a flavour index a sum over flavour will be 



understood, e.g.: 7^ = ^q,7^ and ej^g = X^a^/s- Furthermore, we can use Eq. ( p^ ) 
to express the density asymmetries of the scalars in terms of the ones of the fermions, 
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m ■ 

simplifications are applied also to the BE for IaRs and Y/\p, the following set is obtained: 



and to an excellent approximation we can write 7^" = 7^"- After the same notational 



YAa = (-) I If - 2 ) f + (y^i^ + yAH,.. + yA^) ^-^^^ , (5.7) 




yar, = (^) ( ^ - 2 j 7^ - E (^a.„ + y^n. + yA^) ^-^^^ - yAirf , (5i 

Yar, = [es {z) - 6/ {z)] {^-2\^-y^~,^. (5.9) 




It is possible, and formally straightforward, to add to these equations the appropriate terms 
that allow to extend their validity also in the SE regime, that is for sneutrino masses below 



the bound Eq. (LI). In order to do this, we denote by 7|''^ the set of gaugino-mediated 



reactions with chirality flip on the gaugino line that are responsible for processes that equi- 
librate particle-sparticle chemical potentials.** We also denote by 7^^^ the set of reactions 
induced by the higgsino mixing parameter fifj that enforce the chemical equilibrium con- 
dition Hu + Hd = 0. The thermally averaged rates for these reactions can be written in an 
approximated form as: 

\ = ^, % = ^, (5.10) 

where is the equilibrium number density for one fermionic degree of freedom, while nig 
and fijj in these equations have to be understood as effective mass parameters in which 
all coupling constants as well as reaction multiplicities are reabsorbed. Extension of the 



validity of Eqs. (5.7)- (|5^ to the SE domain is then achieved by adding the following terms 
to the equations for Rb and R^: 

YM, = {YAR,}-yAglf, (5.11) 

Yil = {yar, } - lyAi^ if + \ [y^H^ + y^n) ' (^-i^) 

where the {lAij} above stand for the r.h.s of the corresponding equations (^) and (^). 
Note that since the Rb charge of the term is RB{HuHd) = 2, conserves Rb and 



accordingly there is no term proportional to 7!^^ in Eq. (|5.11 ). Since higgsino equilibration 



involves also the density asymmetry 3^^^ we give below the corresponding C vectors to 



"For M ~ 10* GeV, the soft terms corrections to this approximation 7^/7- — 1 = {A'^ — AB)/M^ can 
be safely neglected. 

**Ref. [14] included a similar term 7mssm in the BE for supersymmetric leptogenesis, corresponding to 
the thermally averaged cross section for the photino mediated process e + e -f^ e + e computed in [45] . 
However, in the total cross section the only contributions that do not vanish in the m-y — >■ limit are those 
that, like e.g. + + br, do not enforce SE. Superequilibrating reactions like + -r^ el + sl 

all vanish in the — >■ limit. 
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express it in terms of the basis of the charge-asymmetries: 

Case I : C^" = — '— (14237, 14237, 14237, 1260, -3915) , (5.13) 
827466 ^ ' ^ ' 

Case II : C^'' = ^ (12469, 16768, 16768, 7056, -21924) . (5.14) 

3 X 162332 

We have of course checked that by increasing the values of rrig and /i^) the results of 
integrating the set of BE given by Eq. (5.7) and Eqs. ( |5.1lD -( p.l2 ) converge to the solutions 



of the usual BE for the SE regime (see Appendix y) . 
5.1 NSE Regime: R-genesis in a simple case 

To highlight the role played by the asymmetries of the two R charges, let us define a simple 
scenario, in which lepton flavour effects play basically no role and thus do not shadow the 
new effects. This scenario is defined by the following two conditions: 

• We assume equal branching fractions for the decays of A*" and of N± into the three 
lepton flavours, that is the defined in Eq. ( p. 14 ) are all equal to | implying = |e 



We assume the regime described in Case I, Section 4.3, in which the Yukawa equi- 
librium condition for the electron holds, and thus the three lepton flavours are all 
treated on equal footing (see the 3x3 upper- left corner in the ^-matrix Eq. ( 4.29| )). 



Given the previous condition, it is then useful to define a 'flavour averaged' lepton 
density asymmetry as: 



a 



With these conditions, the three equations for the flavour charges Eq. ( [5.7D can be 
resummed in closed form into a single equation for B — L: 



(5.17) 



1a,_, = -6 [z) 1^ -2j^ + [yM + y^H^ + y^g ) ^1-^ , (5.16) 

yielding a reduced set of just 3 BE. The 3x3 matrix relating {Ym^ ^Ag, ^A/f,,} *° 
three charge-asymmetries {1ab-l'^a_Rb>^a_Rx} can be readily evaluated from Eq. ( |4.2g| ): 

/ -26348 -6255 4659 ' 

y -11379 12012 -37323 J 
It is now easy to see that in the NSE regime we can rewrite the BE as 

^As-i = SlAiJ^ -iARs> (5-18) 

Yar, = es{z) - 2^ 7jv - (^^a^ + 3^ah. + ^^a^) ^^^^^ - y^-^-f , (5.19) 

Yar^ = [es {z) - ef {z)] " 2^ ^ - 3^A^^ , (5-20) 
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since the difference in the r.h.s. of Eq. ( ^.18| ) gives precisely Eq. ( 5.16| ). Eq. ( 5.1^ ) makes 
apparent how Yar^ and Y/^jig, that in the T — )• hmit keep having non vanishing CP 
asymmetries, are sources of the B — L asymmetry. This result is in fact completely general: 
the only role of the two conditions listed above is simply that of allowing to collapse the 
three equations for Aq, into a single one for Ab-l, while maintaining the BE equations 
in closed form. In particular, it holds also when scattering processes are included (see 
Appendix ^) and independently of the particular (NSE) temperature regime and flavour 
configuration. In short, in the NSE regime the evolution of Ab-l can be always obtained 
from the evolution of 3Ar^ — Ajig, and the final value of ^As-i can be equally well 
obtained from summing the values of the flavour charges asymmetries Yla YAci or from 
the final value of 3Y^f{^ — Ya_r,s- The reason why this happens is simple: by using the 
definitions Eqs. (U)-(^) together with Eqs. (|]2|)-(U) one obtains that 3R^ - Rb = 

— l-B — PQ. Of course, only the PQ fragment of this charge is violated in sneutrinos 
interactions, and from Table |^ we see that this violation is precisely the same than for 
B -L (e.g. for N £Hu we have A{B - L) = -AL = -A{PQ) = -1). Thus, regardless 
of the fact that B — L, Rb and R^ are all independent charges, in the NSE regime the BE 
for 3Ya_Rx ~ YaBb '^ill always coincide with the BE for Y^g_^ = '^Za^^a- 

In our particularly simple case we can take a further step. Let us rewrite the density 
asymmetry y^^g and the combination {y^t + 3^^^ +3^Ag) in the r.h.s of Eqs. ( ^.19|) -( 5^20|) 
in terms of Y/^g_^, ^AR^i ^AJ?^ by means of the A matrix Eq. ( 5.17 ). We can then replace 
Yab-l ^^^Rx ~ ^^Rb and, by using = 7^ we obtain: 

3iAi?,, = [es (z) - ef (z)] " ^ " (9152 FAfi, + 15393 FAiiJ ^ ,(5.21) 

YaRb = 2es{z) - 2^ M - (ll4424yAij, - 24551iyAi?J ^ • (5.22) 

These two equations show that although the asymmetries produced in the two charges 
'iR^ and Rb tend to cancel when taking the difference, their respective washouts are quite 
different, and such a cancellation will never occur. In the general case flavour dynamics 
does not allow to collapse the set of BE to just two equations, but still the same mechanism 
is at work: because of the different washouts, the difference between 3Yr^ and Ypig becomes 
of the same order of these density asymmetries, and so does Yas^l- Perhaps surprisingly, 
we can then expect that by increasing the washouts from a strength of order weak up to 
(not too) large strengths, the flnal value oi B — L will increase. The numerical results in 
the next section confirm this picture. 

In the SE regime instead, things proceed in a different way. Eqs. (5.11)-(5.12) show 
that the BE for Yr^ and Yr^ acquire new washout terms, that are proportional to the 
SE rates, while on the contrary no analogous terms enter the BE for Ya„ Eq. (|5.7| ) or for 
Yab-l Eq. (Hi). Thus, in the SE regime, Eq. ( ^.ISj ) does not hold. One can argue instead 
that, because of the SE washouts, the roles of Ab-l and of 3Aii^ — ARb get reversed, 
since now we have 

3 lAR, - Y^Rb = Yab-, + (yAH^ + yAH,) 7^^ • (5.23) 
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In other words, since SE reactions conserve B — L but violate the R and PQ charges, 
the only source of asymmetry surviving SE is the (thermally induced) Iab-l asymmetry. 
Given that Ai?^ and Ai?^ both contain 'fragments' that carry B number, they do not 
vanish in the SE limit, but are driven to values that are proportional to Ab-l- The 
constants of proportionality are determined by the chemical equilibrium and conservation 
law conditions appropriate for the specific regime and, for example, in Case I are given by 
Yarb = -IYas_l and Yah,, = -^1as-l- 



6. Results 

In this section we quantify the results that are obtained for the baryon asymmetry yield of 
soft leptogenesis when the effective theory described in the previous sections is used, and 
we confront them with what is obtained in the standard scenario, in which SE is assumed to 
hold at all temperatures. Our results are obtained by numerical integration of the BE given 
in Appendix ^ that also include various scattering processes. The comparative results for 
the SE case can be obtained in two formally different, but physically equivalent, ways. A 
first possibility is that of taking the limit m^, //^j — t- oo in the complete BE (given, for 
example, in their basic form in Eqs. ( 5.1l| )- (|5.12| )). A second possibility, that corresponds 



to usual treatments, is to solve only the three equations for the flavour charge-density 
asymmetries with the corresponding A matrix and C vectors obtained under the 
assumption of SE. For the two cases we are analyzing: Case I (/ie,d Yukawa equilibrium) and 
Case II (/ie,(i Yukawa non-equilibrium), we give the corresponding matrices in Appendix^ 
in Eqs. ( |B.13| )-( ^.14| ). Of course, we have verified that both procedures yield the same 
results. 

Some of our results are presented in terms of an efficiency parameter r] defined according 

to: 

Yab-l^z oo) 



V 



2-eY3{z^Q) 



(6.1) 



where e is defined in Eq.( p.l5 ). To single out the new effects that we want to quantify, 
all our results are obtained assuming a configuration of flavour equipartition, with all the 
flavour branching fractions Eq. ( 3.14 ) equal: Pq, = ^, so that flavour effects are basically 



switched off. In all cases, the heavy sneutrino mass is held fixed at M = 10^ GeV, that 
is above the temperature threshold for the validity of the effective theory Eq. ([4.1[). The 



values of the other relevant parameters are: A = 1 TeV, (f)A = ^ and ^ = ^ = 10^^ 
that corresponds to a resonantly enhanced CP asymmetry in mixing. This is obtained for 
2 2? ~ r ~ 2.6 ( o"i'eV ) G^sV. As regards gaugino mass dependent contributions to the CP 
asymmetries from vertex corrections, as was mentioned in Section |2| they are suppressed by 
additional powers of Asusy/M. Given the large value of M that we are using, they remain 
irrelevant even in the cases labeled as the "m^ — )• c« limit" , since in practice nig ~ 10 TeV is 
more than sufficient to enforce SE, and this is the value we are effectively using. Therefore, 
in our regime e is essentially determined only by CP violation in mixing. 

We plot in Figure || the evolution of Ia^.^ with increasing z = M/T. The solid (red) 
lines correspond to the full results obtained in the m-g, Hh ^ OGeV limit, that is when 
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Figure II: Evolution of Yab_l- The solid continuous (red) line depict the complete results in the 
TTig = fif, ^ limit. The dashed (blue) line correspond to the same limit but with all thermal 
corrections to the CP asymmetries neglected. The dotted (black) line gives i^Ai5_i with thermal 
effects when SE is assumed. Panels on the left and right sides are respectively for Case I {he,d 
Yukawa equilibrium) and Case II (/ie,d Yukawa non-equilibrium). Upper and lower panels are 
respectively for rrictf = 0.05 eV and rrictf — 0.20 eV. 



particles-sparticles superequilibrating processes are completely switched off. The dashed 
(blue) lines give the results obtained in the same limit, but when all thermal corrections 
to the CP asymmetries are neglected, and = — = e/2. From all the four panels we 
see that in the NSE regime neglecting thermal corrections is an excellent approximation 
that reproduces with very good accuracy the (sizable) final values of Ias.l- The dotted 
(black) lines give Y^g_^ in the usual treatments which includes thermal corrections and 
also assumes SE, that in our treatment corresponds to taking the limit rrig, //jj — t- oo. 
Panels on the left side refer to Case I discussed in Section 4.3, panels on the right side are 
for Case II discussed in Section iA. We can see that the differences between the situations 
in which the h^^d Yukawa reactions are in equilibrium and when they are out of equilibrium 
are rather mild. Therefore in the following we will concentrate just on results for Case 
I. Upper and lower panels correspond instead to two different strength for the washout 
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processes, parameterized respectively by meg = 0.05 eV and rngfj = 0.20 eV. As it was 
expected from the analysis in the previous section, we see that the stronger the washouts, 
the larger is the gain in efficiency with respect to the SE scenario. 




Figure III: Efficiency factor 77 as a function of the washout parameter nics for Case I {he,d Yukawa 
equilibrium) and different values of rrig ~ fj-a- The red continuous line corresponds to the nig = 
Hh = 100 GeV which is still in the full NSE regime, while the dashed blue line to the same 
limit but with thermal corrections neglected. The red dash-dotted line corresponds respectively to 
■nig — fj,H — 500 GeV, and the black dotted line to SE with nig, Hh — > 00. 



In Figure [D we plot for Case I the efficiency r] defined in Eq. ( |6.lD as a function of 
the washout parameter rries- The red continuous line corresponds to m-g = fj.^ = 100 GeV. 
We have chosen a non-zero value for these parameters because of phenomenological mo- 
tivations, however we have checked that the results are practically indistinguishable from 



those obtained in the m-g = fijj — )• limit and thus, in agreement with Eq. (4.1), the 
evolution still occurs in the full NSE regime. The red dash-dotted line corresponds to 
= fJ-ii = 500 GeV. We can see that in this case SE rates start suppressing the efficiency, 
but are still far from attaining full thermal equilibrium. The black dotted line corresponds 
to the m-g, fijj ^ 00 limit of complete SE. We see that if SE is incorrectly assumed in tem- 
perature ranges where it does not occur, one could vastly underestimate the leptogenesis 
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Figure IV: The final value of Yab-l normalized to the SE result Y£^_^ as a function of nig and 
IJ,H for Case I {h^^d Yukawa equilibrium) and = 0.20 eV. The red continuous line corresponds 
to varying simultaneously both parameters holding nig = fif,. The blue dashed line corresponds to 
varying only rrig in the limit /Ltj} — >■ oo. The green dotted line corresponds to varying only /x^f in the 
limit rrig — ^ oo. 

efficiency. The size of this underestimation is a fast increasing function of the washouts, 
and for particularly large values of mes can reach the two orders of magnitude level. Let 
us also note that for ^ 6 x 10~^ eV, the assumption of SE results in a baryon asym- 
metry of the wrong sign. Graphically, one can see this from the fact that at small values 
of rUeS the black and red lines approximately overlap, and then both change sign around 
TTieff ~ 3 X 10~^ eV. But around rrics ~ 6 x 10~^ eV for the red line there is another sign 
change. This marks the onset of i?-genesis domination; therefore, from this point onward, 
baryogenesis does not proceed through leptogenesis, but rather through i?-genesis. 

In the same figure we have also plotted with the dash blue continuous line the NSE 
results in the approximation of neglecting all thermal corrections to the CP asymmetries. 
By comparing with the full results (red continuous line) we see that for m^s ^ few x 10~^ eV 
thermal corrections give negligible effects. We conclude that in the case of i?-genesis, the 
zero temperature approximation yields quite reliable results. 



-27- 



10 



-10 



10 



-11 



10 



10 



•12 



10 



•13 



10 



-14 



1 1 1 1 1 1 1 

rm,,,=0.2 eV 


1 1 1 1 1 1 


1 1 1 1 1 
B-L ^ 












-. Rb 






Rx ^ 







■ 1- ' 


1 

1 — 
ii' 



10 10^ 

jLi=m's'(Ge\/) 



10 



Figure V: Final values of the charge density asymmetries as a function of nig = /i^j for Case I {h^^d 
Yukawa equilibrium) and rrieff = 0.20 eV. Thick red line: Y/^^ j^; thick blue line: SYaa^ — Y/^Hg] 
thin dashed purple line: Y/^n^ ; thin dotted purple line: Yah • 



In Figure |^ we plot the final value of Y^^ ^^ as a function of different values of nig 
and iJ,fi, normalized for convenience to the value ^ obtained when SE is assumed. 
The results correspond again to Case I discussed in Section 13. In order to enhance the 
impact of the new effects, we have fixed the washout parameter to a rather large value 
TTiefj = 0.20 eV. The red continuous line corresponds to varying simultaneously both SE 
parameters keeping their values equal: rrig = fijj . We see that for m-g = ~ 1 TeV the 
amount oi B — L asymmetry produced by soft leptogenesis can be up to two orders of 
magnitude larger (and of the opposite sign) with respect to what would be obtained in 
the usual approach with SE. SE effects start suppressing the asymmetry around m-g = 

^ 1 TeV. The asymmetry then changes sign around 3 TeV, that marks the transition 
from the i?-genesis to the leptogenesis regime, and eventually around 5 TeV SE reactions 
attain complete thermal equilibrium and Yiig_i^/Y^^ l ^ ^' some interest 

knowing what happens if only one of the two anomalous symmetries U{l)ii or U{1) pq were 
present. While we have not constructed such theories, our BE equations are sufficiently 
general to allow exploring numerically also these cases. The corresponding results are also 
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depicted in Figure The blue dashed line corresponds to the C/(l)R-theory where m-g 
is varied while U{1)pq is broken. The green dotted line corresponds to the alternative 
[/(l)pQ-theory in which rrig — >• oo and only /i^ is varied. From these results we see that 
the real responsible of the large effects is the i?-symmetry, while the effects of the PQ 
symmetry remains qualitatively more at the level of typical spectator effects. A theoretical 
justification of this behavior is not difficult to find, and we will discuss it in the following 
concluding section. 

Some important aspects of the transition from i?-genesis (NSE regime) to leptogenesis 
(SE regime) are highlighted in figure |^ where we plot the final value of the relevant charge 
density-asymmetries as a function of m-g = fJ-a, assuming Case I and nif-s = 0.20 eV. The 
thick solid red line corresponds to Yas-l' while the thin solid blue line corresponds to 
3^A_Rx — YaRb ■ The thin dashed and dotted purple lines display respectively ^aRs and 
^AR^ ■ We see that up to mg = /x^j ~ 100 GeV we have ^As-l — '^^ar^ — ^A_Rs that 
is in agreement with Eq. ( |5.18| ), and thus implies that baryogenesis occurs almost only 
via i?-genesis. As the soft supersymmetry-breaking parameters are increased, SE reactions 
begin to wash out efficiently Y^j^g and ^aHx ^ut the difference 3YAi?x ~ ^a_Rb still remains 
of the order of lAs-i , and i2-genesis still gives the dominant contribution to baryogenesis. 

Around rrig = /Xjj ~ 3TeV all the charge asymmetries change simultaneously their 
sign. This is the benchmark of the onset of the regime in which leptogenesis dominates. 
The only relevant source for generating the density-asymmetries is now the (opposite-sign) 
thermally induced B — L asymmetry, that is not affected by SE washouts, and that is 
feeding (small) asymmetries into all the other charges. In this regime ^A^^ and Y^^^ 
do not have anymore an independent dynamics, and can be simply computed in terms of 
>Afl_z, yielding Y/^r^ = -^Yab-l and Iaa^ = -^^As_i,- 

7. Discussion and conclusions 

The super symmetric seesaw model unavoidably entails the possibility of soft leptogenesis. 
The interest in this possibility relies on the fact that while supersymmetric leptogenesis can 
only proceed within temperature regimes that are in strong tension with the bounds from 
overproduction of gravitinos, typical soft leptogenesis temperatures are sensibly lower, and 
can accordingly relax this tension. However, soft leptogenesis is plagued by the problem 
of a congenital low efficiency, that is related to the cancellation between the asymmetries 
produced in fermions and bosons carrying lepton number. As we have discussed in length, 
this cancellation becomes almost exact in the zero temperature limit. Eventually, finite 
temperature corrections, that break supersymmetry and spoil the cancellation between the 
scalar and fermion CP asymmetries, can rescue soft leptogenesis from a complete failure. 

It should be stressed at this point that the fact that lepton number L commutes with 
supersymmetric transformations, that is that scalar and fermionic members of a supermul- 
tiplet have the same lepton number, plays a crucial role in enforcing the aforementioned 
CP asymmetry cancellation. 

^^Note that since breaks both symmetries, the case of the C/(l)_R-theory is somewhat academic. We 
include it to put in evidence the fundamental role of U(1)r in enhancing the baryon asymmetry. 
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In this paper we have pointed out that in the temperature regime quantified by 
Eq. (0), in which ah reactions that depend on the soft gaugino masses do not occur, 
the early Universe efi'ective theory includes a new i?-symmetry. In soft leptogenesis, this 
i?-symmetry is violated in the out of equilibrium interactions of neutrinos and sneutrinos. 
In particular, i?-number CP asymmetries in heavy sneutrino decays can be defined, and 
constitute important quantities. In fact, given that /^-symmetries do not commute with 
supersymmetry transformations, it is hardly surprising that no cancellation occurs between 
the i?- number CP asymmetries for scalars and fermions. For this reason, a sizable density 
asymmetry for the R charge can develop in the thermal bath, and this asymmetry turns 
out to be the main responsible for the generation of the baryon asymmetry. 

To keep higgsinos sufficiently light, in supersymmetry one needs to assume /i^ ~ T^g, 
and thus when the gaugino masses are set to zero, one must set /i^j — t- as well. In this limit 
the effective theory acquires another quasi-conserved global symmetry, that is a U{1)pq 
symmetry of the Peccei-Quinn type. PQ is also violated in sneutrino interactions and thus 
it also has an associated CP asymmetry. However, since U{1)pq is a bosonic symmetry 
that commutes with supersymmetry, the same cancellation between fermion/boson CP 
asymmetries occurring for lepton number also occurs for PQ. Accordingly, PQ does not 
play an equivalently important role in the generation of the baryon asymmetry. 

In order to make more understandable the previous two remarks, let us start from the 
beginning, by listing the relevant global symmetries of the effective theory. For simplicity 
we concentrate on Case I (/ie,d Yukawa equilibrium). Neglecting lepton flavour, that is 
irrelevant for the present discussion, these symmetries are: L, R, PQ, B and X„c . The flrst 
three L, i?, PQ are violated perturbatively in the interactions of the heavy sneutrinos, and 
all five symmetries are violated by non-perturbative sphaleron processes. In this paper, in 
carrying out our analysis, we have first identified the anomaly free combinations of the five 
charges, that are B — L, Rb and R^, and then we have written down the BE to describe 
their evolution. Here, we want to sketch a different procedure. We first write a set of 
evolution equations for the five anomalous charges, that have the form: 

YAa=SAQ+GAa + GAa- (7-1) 

In this equation S represent the source term for Ya, Q is the (s)neutrino-related washouts 
with all density-asymmetries and signs absorbed, and Q^^ represents the non-perturbative 
EW and/or QCD sphaleron reactions that violate Aq. The latter are reactions of type (i), 
that is fast processes, that eventually will be convenient to eliminate in favor of chemical 
equilibrium conditions. Now, given that B and X„c are good symmetries at the perturbative 
level, they have no CP-violating source term and Sabi ^a^ = (they also do not have 
perturbative washouts, and Qab ~ ^ too). The only source terms thus are ^a^, <Sapq 
and Sah- However, as we already know, in the T — )• limit, for Sa^ have a cancellation 
between the fermion and scalar contributions: + — )• 0. This straightforwardly 
implies that 5^^^ + ^Apq ~^ ^ ^^O' since the sneutrino processes contributing to the CP 
asymmetry for PQ are the same than for L: they are simply multiplied by the appropriate 
PQ charge that is, however, the same for fermion and scalar final states. For the R charge 
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we have instead (x Rf ■ S^^ + Rg ■ S'^^, where Rf^g are respectively the overall R- 
charges of the fermion and boson two particle final state, and thus satisfy Rg = Rf + 2. 
We then straightforwardly obtain that in the T — ^ limit the -R-charge source term does 
not vanish, and is given by Sa^ ^'^A^- Fast in-equilibrium sphaleron processes enforce 
equilibrium conditions between particle densities carrying R charge, and those carrying a 
B and L numbers and, as a result, eventually baryon and lepton asymmetries roughly of 
the same order than the R charge- asymmetry develop. Eventually, with the decreasing of 
the temperature, gaugino mass related reactions will start occurring with in-equilibrium 
rates erasing any asymmetry in the R charge. It is important to notice that when the 
i?-symmetry gets explicitly broken, generalized EW sphalerons reduce to the standard EW 
sphalerons and sphaleron induced multi-fermion operators decouple from gauginos, and 
reduce to their standard B + L violating form. Since gaugino mass reactions as well as all 
other MSSM processes conserve B—L, the asymmetry initially generated through i?-genesis 
will remain unaffected. 

Now that we have identified where the large density asymmetries come from, we can 
complete our procedure by constructing suitable linear combinations of the five equations 
(7.1) for which the sphaleron terms Q^^ cancel out. Since there are only two such terms. 



Gew ^^'^ Gqcd-: '^^^ construct three linear combinations in which only processes of 
type (iii) enter. These are the BE equations for the three anomaly free charges that have 



been discussed at length in Section 4.1. The equilibrium conditions enforced by and 



Gqcd have to be imposed on the system, and to obtain the BE in closed form, the various 
density-asymmetries appearing in the washout terms G must be rotated into the densities 
of the anomaly free charges by means of the appropriate A matrix. 

In this paper, we have not formulated possible alternative effective theories in which 
for example only fi^ = is set to zero, that would correspond to an C/(l)pQ-theory, or 
the alternative case of having just an [7(l)^-theory. However, we have written down a 
set of BE that are sufficiently general to allow exploring numerically the outcome of such 



scenarios. The corresponding results are resumed in Figure |V|, and confirm the crucial role 
played by the R symmetry. In contrast, the effects ascribable to the new PQ symmetry 
arising in the fia ^ limit, that as we have seen are not related with any new large CP 
violating source, remain of the typical size of spectator effects. 

In conclusion, supersymmetry offers different ways to explain the cosmic matter- 
antimatter asymmetry. The asymmetry could be directly generated in baryon number 
since, although severely constrained, EW baryogenesis has not been ruled out yet. Alter- 
natively, the asymmetry could be initially generated in lepton number, through supersym- 
metric leptogenesis [13] or through soft-leptogenesis if it occurs below T ~ lO'^ GeV. The 
main finding of our paper is that there is also a third, previously unnoticed, possibility. 
That is that the asymmetry can be first generated in the new R charge that appears in the 
effective theory for supersymmetry when the Universe temperature is above T ~ 10^ GeV, 
and then transferred to baryons via generalized EW sphalerons. 



We are concentrating here on the role and fate of the 7?-symmetry. However, given that eventually also 
the PQ symmetry gets explicitly broken, higgsinos decouple from sphalerons as well. 
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A. Thermal factors 



In terms of the dimensionless evolution parameter z = M/T the thermal factors appearing 
in the expressions of the decay CP asymmetries Eqs. ( |3.9| )-( |3?T^ ) read: 



c'^^iz) 



c^{z) + cf{z)' 

where, in the approximation in which N± decay at rest, 

C^(z) = {1-Xi- XgJX{l,Xi,XgJ [l - fP] 



1 - fJ 



with 



C^{z) = \{l,XH^,Xj)[l + fl] 
Xa[.z) - 



eq 



M2 ' 

A(l,x,y) = A/(l + x-y)2 -4x. 

The Bose-Einstein (s) and Fermi-Dirac (/) equihbrium distributions are: 

1 



eq 



f 



eq 



1 



where 
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Finally, the thermal masses for the relevant scalar and fermion particle species are [15]: 

(A.IO 



xj = 2xi 



1 /3 



:92 + IdY ) : 



(A.ll 



where g2 and gy are the SU{2) and U{1) gauge couplings, and is the top Yukawa 
coupling, renormalized at the appropriate energy scale. 
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B. Boltzmann Equations 



In this Appendix we present the Boltzmann equations that must be used for numerical 
studies of soft leptogenesis when the heavy sneutrino masses satisfy the condition Eq. . 
We also include the SE reactions and 7^^^ defined in Eq. ( 5.1C| ) , that extend the validity 
of our BE to all temperatures. 

The Boltzmann equations which describe the evolution of RH neutrino and sneutrino 
densities are: 



Yn = -[^-i) (7. + 47f ) + 47« + 47f ) + 27^^ + 47^^) , (B.l) 
= - - 2) (f + 37.. + 27^^ + 27f ) + 27r) + 7f ^ + 27^^) , (B.2) 



where the time derivative is defined as 1" = sHz^, s is the entropy density, and H = H{z) 
is the Hubble parameter. We have defined = + Y^ , and the reaction rates 7 
without a flavour index a are always understood to be summed over all flavours. For the 
evolution of the flavour charges Ia^ we have 



Ya^ = -[Ea + Eo^] , (B.3) 



where 



- (7^^"^ + 27^^" + 27f + 27^^" + 7f 3^A.. 

- U'"- + 7f + + 7f + 7f + U''^^) y^-^ 

\ ^ N+J 

I (5)fc , {7)k , 1 (6)fc Y^ \ / 

- hi + 7i + ^Yt I r^AH„ - ^A//„ 

+lf{yM-yM^, (B.4) 
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and 



N 



s.a 

'n 



2 I 



1 Y~ 

N 



N+ / 



+ 27^13^ 
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(9)a -^AT 
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2y 
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-if (y^j^ - y^ 



(B.5) 



The 3^A appearing in these equations are defined in Eq. ( |5.4| ), while the SE reaction rate 
if has been defined in Eq. ( ^.101 ). For the decay reaction densities we have: 



7~ 

'n 



7~ 

' AT 



f,a . s,a 

7^ +7^ , 



AB 
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(B.6) 



where A and i? are taken to be real. For values M ~ 10^ GeV the higher order terms in 
the soft parameters can be safely neglected. 
The scattering processes considered are 



Reaction 
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where for convenience we have listed the corresponding changes of the R-charges in each 
process. The reduced cross sections for the processes listed above can be found in ref. [14]. 
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The BE above do not include the CP asymmetries of top and stop scatterings. Strictly 
speaking, when scatterings are included, for consistency one should include also the corre- 
sponding CP asymmetries. However, in soft leptogenesis this cannot be done in a straight- 
forward way because thermal factors for the scattering CP asymmetries constitute a new 
set of non trivial quantities. Fortunately, in the strong washout regime for leptogenesis, the 
effects of CP asymmetries in scattering have been found to be subdominant with respect 
to CP asymmetries in decays [46], and since in this paper we focus precisely on strong 
washouts, neglecting the scattering CP asymmetries is justified. 



The BE for the evolution of Rb and Ry., defined in Eqs. (4.6) 



are: 



7f 



eff 



' yAa + ^[y^H+yAH.r 



(B.7) 
(B.8) 



where again the SE rates 7- and 7^^^ have been also included. Fa and Fa are given by: 
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For Ga and Ga we have: 



Go 
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As we have explained, 




yAH^- 



(B.12) 



with the inclusion of 7- 



and 7^'': our BE are valid at all 



temperatures. To verify this, we have compared the results obtained with the complete 
BE given above, with what is obtained by integrating the set of BE specific for the SE 
regime, that reduce to the equations for the neutrino and sneutrino abundances Eq. ( |B.1| ) 
and Eq. ( |B.2D plus the three equations for the flavour charges Eq. ( |B.3| ). Of course, one 
also has to use the matrices and C^" vectors appropriate for the SE limits of the two 
cases that we have been studying (recalling also that = 2A^ and G^^ = 2C^"). For 
Case I of Section 4.S we have: 

/ 



A' 
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9 X 237 
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C 



237 " 



(B.13) 



that, incidentally, coincides with the matrix given in [13] for the case of all Yukawa couplings 
in equilibrium. The matrix for Case II of Section 4.4 is given in [13], and is rewritten below 
for convenience: 
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(37, 52, 52) . (B.14) 
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